Three-dimensional topological insulators in the octahedron-decorated cubic lattice 
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We investigate a tight-binding model of the octahedron-decorated cubic lattice with spin-orbit 
coupling. We calculate the band structure of the lattice and evaluate the Z2 topological indices. 
According to the Z2 topological indices and the band structure, we present the phase diagrams of 
the lattice with different filling fractions. We find that the (1; 111) and (1;000) strong topological 
insulators occur in some range of parameters at 1/6, 1/2 and 2/3 filling fractions. Additionally, the 
(0; 111) weak topological insulator is found at 1/6 and 2/3 filing fractions. We analyze and discuss 
the characteristics of these topological insulators and their surfaces states. 

PACS numbers: 73.43.-f, 71.10.Fd, 73.20.-r, 72.25.-b 



I. INTRODUCTION 

Usually, different phases of matter can be classified 
using Landau's approach according to their underling 
symmetries pj. In 1980s, the discovery of the quantum 
Hall effect changed physicists' viewpoint on the classifi- 
cation of matter [2|. The quantum Hall states can be clas- 
sified by a topological invariant, now named the TKNN 
number [3I (equivalent to the first Chern number), which 
is directly connected to the quantized Hall conductiv- 
ity, but they have the same symmetry. Since the Hall 
conductivity is odd under time reversal, the topological 
non-trivial quantum Hall states can only occur when time 
reversal symmetry is broken, which is performed by a 
magnetic field. In 1988, Haldane also proposed a time 
reversal symmetry broken toy model without a magnetic 
field to realize quantum Hall states [4]. All the quantum 
Hall states have a gapped band structure in bulk and chi- 
ral gapless edge states that are topologically protected. 

Recently, the promising prospect of spintronics in tech- 
nology stimulates physicists to generate spin current. 
Quantum spin Hall effect was proposed to create spin 
current [il, |6[. The quantum spin Hall states are non- 
trivial topological phases with time reversal symmetry, 
which have a bulk gap and topologically protected gap- 
less helical edge states. For the above reason, the quan- 
tum spin Hall states also called topological insulators. 
Two-dimensional topological insulators are characterized 
by a Z2 topological index u = 0, l[7|. For a non-trivial 
topological insulator the topological index has a value 
V = 1 while = for a trivial band insulator. There- 
fore, a topological insulator always has a metallic bound- 
ary when placed next to a vacuum or an ordinary band 
insulator because topological invariants cannot change 
as long as a material remains insulating. The remark- 
able metallic boundaries of topological insulators may 
result in new spintronic or magnetoelectric devices and 
a new architecture for topological quantum bits. In 
quantum spin Hall phases, the spin-orbit coupling plays 
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the role of the spin-dependent effective magnetic field. 
The first real material, a HgTe quantum well, support- 
ing two-dimensional topological insulators was predicted 
by Bernevig, et a/. [8] and experimentally conformed by 
Konig et a/. [9.]. 




FIG. 1: (Color online), (a) The octahedron-decorated cubic 
lattice which can be obtained by replacing every lattice site 
of a cubic lattice with an octahedral cluster as shown in (b). 
(c) The three-dimensional Brillouin zone and high symmetry 
points, (d) The two-dimensional Brillouin zone of a slab with 
two 001 surfaces. 

Soon after the quantum spin Hall insulator was dis- 
covered, time-reversal invariant topological insulators 
were generalized to three dimensions p^Ql - [l2| . Three- 
dimensional time-reversal invariant band insulators are 
classified according to four Z2 topological indices 
(vq] ^1^21^2) with Vi = 0, i[1Q|. In three dimensions, 
the time-reversal invariant band insulators can be clas- 
sified into 16 phases according to the four Z2 topolog- 
ical indices. A band insulator with = 1 is called a 
strong topological insulator (STI), a band insulator with 
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z/Q = and at least one non-zero = 1, 2, 3) is called a 
weak topological insulator (WTI), while an ordinary triv- 
ial band insulator has an index (0;000). For an STI 
phase, the surface states have an odd number of Dirac 
points, which are topologically protected and for a WTI 
or trivial band insulator phase, the surface states have 
an even number of Dirac points. Fu and Kane firstly 
predicted that Bii-x Shx supports a three-dimensional 
topological insulator [13], which was conformed experi- 
mentally by Hsieh, et al. in 2008 [l^- Later, Bi2Se3 
was discovered to be a three-dimensional insulator exper- 
imentally as a second generation material [l5| , which also 
was supported by theoretical calculations [TBI [HI- Ad- 
ditionally, reference [HI also predicted that Bi2Te3 and 
Sb2Te3 are second generation materials supporting three- 
dimensional topological insulators. The later experimen- 
tal studies on Bi2Te3[i3-[i3 and Sb2Te3il^] identified 
their topological band structures. 

To help experimental physicists find more topologi- 
cal insulator materials, theoretical physicists have in- 
vestigated several models that support non-trivial topo- 
logical insulators. Theoretical studies have demon- 
strated that, within the tight-binding approximation 
and with the spin-orbit coupling, the honeycombra, 
kagome^20i]. checkerboard [2 li] . decorated honeycomb ^22|, 
Lieb[23[, and square-octagon [2^ lattices support two- 
dimensional topological insulators and the diamond [Tq|, 
pvrochlore[25|. and perovskite[23] lattices support three- 
dimensional topological insulators. 

In this paper, we shall show that a new lattice, the 
octahedron-decorated cubic lattice as shown in Fig.l (a), 
supports three-dimensional topological insulators with 
the spin-orbit coupling existing. This lattice can be 
regarded as a three-dimensional generalization of the 
square-octagon lattice [22]. We find that this model sup- 
ports STI and WTI phases for 1/6 and 2/3 filling and STI 
phases for 1/2 filling as well as ordinary band insulator 
and metal phases. 



II. MODEL 



We consider the octahedron-decorated cubic lattice as 
shown in Fig(T] (a), which can be obtained by replacing 
every lattice site of a cubic lattice with an octahedral 
cluster as shown in Fig(T]^b). This lattice has a unit cell 
with six different lattice sites as denoted in Figdjb) so 
that it contains six sublattices. Here, we assume that the 
distance between the centers of two nearest-neighbor oc- 
tahedral clusters is a, which is the same with the lattice 
constant of all sublattices, the distance of every lattice 
site of an octahedral cluster from its center is a/4, and 
the distance of two nearest-neighbor lattice sites in dif- 
ferent octahedral clusters is a/2. With the tight-binding 
approximation, we can write the second quantized Hamil- 



tonian of the lattice as follows. 
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(1) 



where Qct is the annihilation operator destructing an 
electron with spin cr on the site of the octahedron- 
decorated cubic lattice, (i, j) represents nearest-neighbor 
hopping in the same octahedral cluster with amplitude t 
and [i, j] denotes nearest-neighbor hopping between two 
different octahedral clusters with amplitude ti. 
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FIG. 2: (Color online). Phase diagrams of the octahedron- 
decorated cubic lattice for (a) 1/6 filling, (b) 1/2 filling, and 
(c) 2/3 filling. Here, BI denotes a trivial band insulator; STIi 
and STI2 denote (1; 111) and (1;000) strong topological insu- 
lators, respectively; WTI denotes a (0; 111) weak topological 
insulator; and M denotes a metal phase. 

In momentum space, the Hamiltonian ([1]) can be 
represented by Hq = Eka^L'^k^^ka with ^ka = 
(cika, C2ka, C3ka, C4ka, Cska, Ceka) , which are Ordered ac- 
cording to the sequence denoted in FigH^b). Here, H^"^ 
takes the following form, 
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Since Hq is spin-decouphng, is spin-independent, i.e. 
it is the same for both spin- up and spin-down electrons. 
FiglUc) shows the first Brillouin zone of the octahedron- 
decorated cubic lattice. The spectrum of Eq.([2]) with 
ti = t is calculated and shown in Figl3l^d). The spec- 
trum contains six bands which come from the six sites 
in every unit cell. A gap exists between the first and 
second bands. The second, third and fourth bands touch 
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together at points T^R and M. The third, fourth and 
fifth band touch at point X, near which a Dirac cone 
occurs. Five bands including the second, third, fourth, 
fifth and sixth bands meet at point F. Along the T ^ R 
line in momentum space, the second and third bands are 
degenerate and the fifth and sixth bands are degenerate. 

Now, in order to find non-trivial topological insulators 
in the octahedron-decorated cubic lattice, we proceed 
to introduce the spin-orbit interactions between next- 
nearest-neighbor sites as follows, 

where {{i^j)) represents two next-nearest-neighbor sites 
i, j, and Aso is the amplitude of spin-orbit coupling of 
the two next-nearest-neighbor sites, a = (cr^^, cr^, cr^) 
is the vector of Pauh spin matrices, djj^ are the two 
nearest neighbor bond vectors traversed between sites i 
and j with S\djj x d|^|/a^ = 1. In momentum space, 
the Hamiltonian for spin-orbit coupling (|3|) can be ex- 
pressed as Hso = Ek^k'^k^^k with ^k = (cikt,C2kt, 

C3kt7 C4kt7 C5kt7 C6kt7 Clki, C2ki, Csk^, C4k|, C5k|, Cek^)^. 

Since does not decouple for the two spin projec- 

tions, it is a 12 X 12 matrix. In momentum space, the 
total single particle Hamiltonian is Hk = + ^k^- 
The bands and eigenstates can be obtained by exactly 
diagonalizing Hk- 



III. THREE-DIMENSIONAL TOPOLOGICAL 
INSULATORS 

The classification of three-dimensional topologi- 
cal insulators is presented in Ref . p^of . For three- 
dimensional lattices there eight distinct time reversal 
invariant momenta (TRIM), which can be expressed 
in terms of primitive reciprocal lattice vectors as 
ri=(ni,n2,n3) = (^1^1 + n2b2 + n3b3)/2 with rij = 0, 1. 
Three-dimensional topological insulators can be distin- 
guished by four Z2 topological invariants (^^o; ^1^2^3)7 
which are defined as (—1)^° = nnj=o 1 ^^i^2n3 and 

( — 1)^*-^'^'^ = rinj^i=0,l;ni = l ^^i^2n3 7 whcrC Snin2n3 = 

^det[w(rn^ri2ns)]/^^[^{^nin2ns)] = ±1- Here the 
unitary matrix w is defined as Wij(k) = {ui{—'k)\Q\uj{'k) 
with O being the time reversal operator and |iij(k)) 
being the Bloch wave functions for occupied bands. Fu 
and Kane have found a simple method to identify the Z2 
invariants for the system with the presence of inversion 
symmetry [l3|. In this case, Sn^n2n3 can be calculated by 
^nmsna = nm=i (Fnmsna ) , whcrc TV is the number 
of occupied bands and ^2m(Fnin2n3) = ±1 is the parity 
eigenvalue of the 277ith occupied band at ^nin2n3' 

Our 

model is inversion symmetric so we will adopt this 
method to evaluate the Z2 invariants = 0,1,2,3). 
We select the center of an octahedron in the lattice 
as the center of inversion, then the parity operator 
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FIG. 3: (Color online). Band structures of the octahedron- 
decorated cubic lattice for various parameters ti and Aso. 
Here, the horizontal axis represents the wave vectors along the 
path in the first Brillouin zone indicated by the red lines in 
FigH^c). (a) ti = t, Aso = 0.5t, (b) ti = -t, Aso = -0.5t, (c) 
ti = t, Aso = t, (d) ti = t, Aso = 0, (e) ti = 3t, Aso = -0.4t, 
(f) ti = -3t, Aso = 0.4t, (g) ti = 3.2t, Aso = -0.2t, (h) ti = 
t, Aso = 0.2t, (i) ti = t, Aso = -0.2t, (j) ti = -t, Aso = 0.2t, 
(k) ti = 2.5t, Aso = -0.2t, and (1) ti = 0.5t, Aso = 0. 



acts as 7^[V^i(r), V^2(r), V^3(r), V^4(r), V^5(r), V^6(r)]^ = 
[V^6(-r), V^4(-r), V^5(-r), V^2(-r), V^3(-r), V^i(-r)]^, 
where [V^i (r), V^2(r), V^aW, V^4(r), V^sW, V^eW]^ is the 
six-component wave function. Taking Fourier trans- 
formation, we can write the six-component wave 
function as [V^i(r), V^2(r), V^aW, V^4(r), V^sW, V^eW] = 
Ek [^1 (k) , ^2 (k) , (/)3 (k) , (/)4 (k) , (/>5 (k) , (k)] e^^- _ and 
the parity operator as 7^ = Xlk '^■^^e 
Then, in momentum space, we obtain the equa- 
tion Pk[<^i(k),<^2(k),<^3(k),(A4(k),05(k),06(k)]^ = 

[M-^), M-^), M-^), M-^), M-^), M-^)f- 

Considering the degree of spin, we can express the parity 
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operator at the time reversal invariant momenta r^^^2n3 
as follows, 
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where the 4x4 matrix is the unit matrix in spin space. 

We diagonalize the total single-particle Hamiltonian 
and calculate the Z2 topological invariants for dif- 
ferent filling fractions. We find that no n- trivial topo- 
logical insulators exist for 1/6,1/2 and 2/3 filling while 
only metal phase occurs for 1/3 and 5/6 filling. Thus, 
we will focus on and discuss the cases with 1/6, 1/2 and 
2/3 filling fractions in the following part of the paper. 
We identify phases for different parameters ti and Aso 
with 1/6, 1/2 and 2/3 filling fractions and draw phase di- 
agrams as shown in Fig 121 Figsl2fa), [2fb) and[2fc) show 
the phase diagrams for 1/6,1/2 and 2/3 filling, respec- 
tively. For 1/6 and 2/3 filling, there are (1; 111) and 
(1; 000) STI phases, (0; 111) WTI phase as weh as trivial 
band insulator and metal phases. For 1/2 filling, there 
are (1; 111) and (1;000) STI phases, trivial band insula- 
tor and metal phases except (0; 111) WTI phase. 

To clearly manifest the bulk band structure of different 
phases for various filling factions, we calculate the bulk 
energy bands for several cases with different parameters 
ti and Aso, which are shown in Fig 131 In order to in- 
vestigate the characteristics of surface states for various 
phases, we evaluate the energy bands in a slab geom- 
etry with two 001 surfaces. The Brillouin zone of the 
slab is shown in Figdjd). The energy bands are present 
along lines that connect the four surface TRIM as shown 
in FigHl With the assistance of the bulk energy bands 
shown in Fig Inland the two-dimensional energy bands for 
a slab shown in FigHl we will sequentially analyze various 
phases, identify three-dimensional topological insulators, 
and discuss their characteristics for 1/6, 1/2 and 2/3 fill- 
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FIG. 4: (Color online). Band structures of a slab with two 
001 surfaces for various parameters ti and Aso- Here, the 
horizontal axis represents the wave vectors along the path 
in the surface Brillouin zone indicated by the red lines in 
FiglHd). (a) ti = t, Aso = 0.5t, (b) ti = -t, Aso = -0.5t, (c) 
ti = t, Aso = t, (d) ti = t, Aso = 0, (e) h = 3t, Aso = -0.4t, 
(f) ti = -3t, Aso = 0.4t, (g) ti = 3.2t, Aso = -0.2t, (h) ti = 
t, Aso = 0.2t, (i) ti = t, Aso = -0.2t, (j) ti = -t, Aso = 0.2t, 
(k) ti = 2.5t, Aso = -0.2t, and (1) ti = 0.5t, Aso = 0. 



FiglJl^a) shows the phase diagram of the octahedron- 
decorated cubic lattice for 1/6 filling. In this case, the 
(1; 111) and (1; 000) STI phases are discovered. The non- 
trivial STI phases have a gap between the first and sec- 
ond bands as shown in Figl3fa) and[3fb) corresponding 
to (1; 111) and (1; 000) STI phases, respectively. We note 
that for 1/6 filling there is only one Dirac point on TRIM 
as shown in FiglH^a) and (b), that is, only a pair of robust 
spin- filtered states exists. We also find a (0; 111) WTI 
phase for 1/6 filling e.g., as shown in Figl3l^c). Fig|4l^c) 
shows the surface states for a (0; 111) WTI phase that 
has two Dirac points between the first and second bands 
on TRIM. We note that trivial band insulators occur for 



smaller ti and smaller Aso parameters, which is easily 
understood for when ti and Aso approaches to zero the 
lattice becomes separated octahedral clusters. For a triv- 
ial band insulator there is a gap between the first and sec- 
ond bands as shown in Figl3l(d), but there are not surface 
states as shown in Fig|4l^d). For a metal phase, the gap 
vanishes. 
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B. 1/2 filling 

FigEfb) shows the phase diagram of the octahedron- 
decorated cubic lattice for 1/2 filling. For 1/2 filling, 
(1; 111) and (1;000) STI phases, trivial band insulators, 
and metal phases occur, but WTI phases are not found. 
Figsini^e) and [3ff) show the band structure for (1; 111) 
and (0; 111) STI phases, respectively. We can find from 
these diagrams that a gap opens between the third and 
fourth bands. For STI phases, there only one Dirac point 
on TRIM as shown in FigHJ^e) and (f). For trivial band 
insulators, there is also a gap between the third and 
fourth bands as shown in Figl3l^g), but even number of 
Dirac points exist on TRIM as shown in Fig|4l^g). For 
smaller ti and smaller Aso, a metal phase occurs except 
a special point ti = and Aso = 0- For ti = and 
Aso = 0, the second, third and fourth bands are degen- 
erate and become a flat band, which means that electrons 
are localized. In other works, the system with ti and Aso 
for 1/2 filling is a trivial band insulator. However, a tiny 
change from ti = and Aso = for parameters ti and 
Aso makes the flat band become three dispersive bands 
that are crossover each other, then the lattice with three 
bands occupied becomes a metal. 

C. 2/3 filling 

Figl2fc) shows the phase diagram of the octahedron- 
decorated cubic lattice for 2/3 filling. We note that, simi- 
lar to 1/6 filling, (1; 111) and (1; 000) STI phases, (0; 111) 
WTI phase, trivial band insulator, and metal phase oc- 
cur in different ranges of parameters ti and Aso- Fig 131(1) 
shows the band structure for ti = Aso = —0.2t at 
which a (1; 111) STI phase occurs. We can find that a 
gap opens between the fourth and fifth bands as shown in 
FiglSKi). There is an odd number of surface states which 
traverse the gap as shown in FigHJ^i). For the (1;000) 
STI phase, the similar characteristics are exemplified in 
FigsEJj) andlH^j). The (0; 111) WTI phase is found as 



weh. FiglSl^k) and Figgl^k) show the (0; 111) WTI phase 
has a gap between the fourth and fifth bands and an even 
number of surface states traversing the gap. For smaller 
ti and smaller Aso, the system for 2/3 filling is a trivial 
band insulator, which is feathered by a gap between the 
fourth and fifth bands combined with an even number of 
surface states traversing the gap as shown in Figl3](l) and 
FigHJl), respectively. 



IV. CONCLUSION 

In summary, we have shown that the octahedron- 
decorated cubic lattice with spin-orbit coupling supports 
three-dimensional topological insulators at 1/6,1/2 and 
2/3 filling fractions. For 1/6 and 2/3 filling, (1; 111) and 
(1; 000) STI phases, (0; 111) WTI phase, trivial band in- 
sulator, and metal phase are found, while for 1/2 fill- 
ing, (1; 111) and (1;000) STI phases, trivial band insu- 
lator, and metal phase occur except (0; 111) WTI phase. 
We have calculated the band structure and surface band 
structure for the tight-binding model of the octahedron- 
decorated cubic lattice with spin-orbit coupling and eval- 
uated the Z2 topological invariants. We have analyzed 
and discussed the characters of the band structures and 
the surface states of different phases. Although the 
octahedron-decorated cubic lattice we considered is a toy 
model, our study points out an alternative path to search 
for real topological materials. On the other hand, it 
might as well be built from optical lattices due to their 
diversity and controllability. 
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